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INTRODUCTION 


If  energy  degradation  (Coulomb  loss)  and  pitch-angle  diffusion  can  be 
neglected,  then  the  radial  transport  of  magnetospheric  ring-current  ions  is 
governed  by  an  equation  of  the  form 

3f/3t  -  L2(3/3L)[(Dll/L2)0?/3L)]  -  (f /r  ) 

-  (3/34)[DM(3f/94)l  -  (f/rq),  (1) 

where  f  is  the  phase-space  density  at  fixed  M  and  J  (first  two  adiabatic  in¬ 
variants),  D^l  is  the  diffusion  coefficient  for  transport  in  L  (dimensionless 
shell  parameter) ,  $  is  the  third  adiabatic  invariant  (inversely  proportional 
to  L) ,  and  is  the  ionic  lifetime  against  charge  exchange.  Jentsch  [1984] 
has  described  a  method  for  obtaining  approximate  steady-state  solutions  of 
this  equation  when  is  exactly  proportional  to  a  fixed  power  (p)  of  L, 
i.e.,  when  is  exactly  proportional  to  a  fixed  power  (4-P)  of  <$.  The  pur¬ 
pose  of  the  present  work  is  to  describe  an  alternative  method  for  obtaining 
time-dependent  as  well  as  steady-state  solutions  of  (1)  while  permitting  the 
dependence  of  upon  L  (or  equivalently,  of  upon  $)  to  deviate  somewhat 
from  a  strict  power  law. 


The  alternative  method  is  highly  advantageous  because  even  the  simplest 
dynamical  models  for  magnetospheric  radial  diffusion  lead  to  diffusion  co¬ 
efficients  Dtt  that  deviate  in  fact  from  strict  power  laws  except  in  certain 
limits.  For  example,  the  standard  model  [e.g.,  Cornwall ,  1972;  Schulz ,  1983] 
for  charged-particle  diffusion  in  a  dipolar  magnetic  field  leads  to  a  diffu¬ 
sion  coefficient  of  the  form 


1 X 10" 10L10  day"1 


LL  (M/y2VZM0)2[2D(y)/T(y)]2  +  10"6L4 

+  7x  10"9(Q(y)/180D(y)]2L10  day"1 


(2) 


where  y  is  the  sine  of  the  equatorial  pitch  angle  y  is  the  ratio  of  rela¬ 
tivistic  mass  m  to  rest  mass  m^,  Z  is  the  integer  that  specifies  charge 
state,  and  Mq  a  1  GeV/gauss.  The  auxiliary  functions  Q(y),  D(y),  and  T(y)  in 


(2)  are 

well  approximated  [Schulz  and  Lanzerotti,  1974,  pp. 

20,  21,  44;  Da- 

vidson. 

1976] 

by 

the  algebraic  expressions 

Q(y> 

a 

-  27.12667  -  45.39913y4  +  5.88256y8, 

(3a) 

D(y) 

=s 

0.4600577  +  0.1Q66154y3/4  -  0.1997662y, 

(3b) 

and 

T(y) 

s 

1.3801730  -  0.6396925y3/4. 

(3c) 

A  further  complication  is  that  y  and  y  in  (2)  typically  vary  with  L  at  fixed 
M  and  J.  The  variation  of  y  is  given  approximately  [Chen  and  Stern,  1975]  by 

y“2  *  1  +  1.38048X  -  0.030425X4/3  +  0.10066X573  +  [X/2T(0)]2,  (4) 

1/2 

where  X  =  (La/8mQ^iM)  J,  £  is  the  radius  of  the  earth,  and  p.  is  the  earth's 
magnetic  moment.  The  variation  of  y  with  L  is  given  by 

Y2  -  1  +  (2tiM/L3a3y2m0c2),  (5) 

where  c  is  the  speed  of  light.  The  limiting  cases  X  »  0  (J  *  0)  and  X  *  = 

2 

(M  *  0)  correspond  to  y  *  1  and  y  *  0,  respectively,  but  X  *  »  implies  M/y 
2  2 

■  J  La/32fAmg[T(0) ]  upon  evaluation  of  the  indeterminate  form. 


TRANSFORMATIONS 

The  factors  [ 2D(y) /y2T(y) ] 2  and  [Q(y)/180D(y) ] 2  in  (2)  vary  approximate¬ 
ly  as  powers  of  L  [Schulz  and  Lanzerotti,  1974,  pp.  91,  93].  The  exponents 
of  L  are  given  by 


(6a) 


a  r 
_ 1 1  n 
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(6b) 

respectively,  where  [Schulz  and  Lanzerotti.  1974,  pp.  20,  21] 


/  3  InjA  ,  _  Y(y), 

V d  ln  wM,j  4T(y} 


0,  y  -  1 
-1/2,  y  -  0 


(7a) 


and 


Y(y) 


(y')~2  T(y')  dy' 


(  0,  y  *  1 

6[T(y)  -  2D(y) ]  -  {  (7b) 

(  2T(0),  y  -  0 


The  exponents  of  L  implied  by  (6)  remain  approximately  (but  only  approximate¬ 
ly)  independent  of  L  for  fixed  (but  nonvanishing)  M  and  J.  They  remain 
strictly  independent  of  L  only  for  y  »  1  and  for  y  *  0,  and  even  in  these 
cases  the  form  of  D  specified  by  (2)  is  not  strictly  proportional  to  a 

LL 

fixed  power  of  L,  e.g.,  proportional  to  with  (3  a  function  of  M  and  J  only. 

The  transformation  z  ■  (P  -  3)ln  L  proposed  by  Jentsch  [1984]  brings  (1) 
into  an  equation  of  the  form 


3w/3t  -  (p  -  3)2L'2Dll[(32/3z2)  -  (l/4)]w  -  (w/-q) , 

where  w  *  ^(P-3)/2|,  if  D  x  LP  for  fixe<i  p,  Xhe  transformation 


(8) 


/; 


a  >2  °L'[  iL 


L  (d«/dL) 


r # 

X  « 


d<S 


(9) 


introduced  in  the  present  work  brings  (1)  into  an  equation  of  the  form 


3f/3t  -  (L4/DLL)(32f/3;2)  -  (f / t  ) 


(10) 


without  recourse  to  the  assumption  that  oc  for  fixed  {3.  The  factor 

0  / 

l  (d$/dL)  in  (9)  is  a  constant,  since  4  oc.l/L,  and  the  factor  L  /D  in  (10) 

LL 

is  proportional  to  for  the  same  reason.  The  fixed  limit  of  integration 

L  ■  «>  (S'  ■  0)  in  (9)  lies  outside  the  domain  of  validity  of  (1)  and  thus 
requires  the  integrand  of  (9)  to  be  evaluated  by  analytical  extrapolation 
for  L  >  L^.  However,  it  follows  from  (9)  that  (,  *  ((3-3)  if 

Dll  *  L*3  for  fixed  (3  >  3,  and  no  other  fixed  limit  of  integration  in  (9) 
would  lead  to  such  a  simple  result  for  The  form  of  (10)  suggests  a  time- 
dependent  solution 

CO 

f(L,t)  -  Ij L)  +  ^  an(t)gn(L)  •  (11) 

n»0 

in  which  f  (L)  is  the  steady-state  solution  of  (10)  and  the  g  (L)  are  the 
x  n 

4  2  2 

eigenfunctions  of  the  operator  A  a  -  (L  /DTT)(3  /9(,  )  +  (1/t  ),  corresponding 

Lia  q 

(respectively)  to  the  eigenvalues  The  expansion  coefficients  an(t)  are 

thus  given  by  an(c)  *  an(0)exp(-\nt)  if  the  transport  coefficients  and  bound¬ 
ary  conditions  are  time-independent  for  t  >  0. 

STEADY  STATE 

The  steady-state  solution  f  (L)  thus  satisfies,  according  to  (10),  the 

CO 

equation 

(d2fa3/d42)  -  (DLL/L4rq)f  -  0,  (12) 

subject  to  the  inner  boundary  condition  that  f^d^)  =■  0  at  the  top  of  the 
atmosphere  and  the  outer  boundary  condition  that  f„(L^)  correspond  to  the 
phase-space  density  at  the  inner  edge  of  the  plasma  sheet.  The  exact  solu¬ 
tion  of  (12)  is  expressible  (cf.  Schiff,  1955  ,  p.  187;  Walt ,  1970,  p.  ■*!•*] 


in  terms  of  modified  Bessel  functions  of  fractional  order  if  Dtt/L  t  is 

LL  cj 

exactly  proportional  to  a  fixed  power  (p)  of  4*  This  fact  suggests  a  modi¬ 
fied  WKB  approximation  of  the  form 


f  (L)  *  f  (L) 

a  a 


's1/2i<vdu.)1M  ys>v8o>  -  vs>iA>  , 

sl/2vyDu.>!/4  -  vAvv  *  1 


for  the  steady-state  solution  of  (1)  and  (10),  where 


3  f  (DL  L  )1/2ta  )4rq(L  }r1/2 


[Tq(L)DLL  ]_1/2  dL' 


r* 

I 

J  o  q 


>W1/2  d*'.  (14) 


The  name  "modified  WKB  approximation"  is  suggested  by  the  appearance  of  modi¬ 
fied  (rather  than  ordinary)  Bessel  functions  in  (13)  as  a  consequence  of  the 
negative  (minus)  sign  in  (12).  The  optimal  order  v  of  the  modified  Bessel 

A  A  _  _ 

functions  1^(0)  and  Ky(0)  in  (13)  is  given  by  v  »  l/(p+2),  where  p  is  a  rep¬ 
resentative  value  of 


p  =  (d  In  D  /d  In  4)  “  4(d  In  L/d  In  4)  -  (d  In  t  /d  In  4) 

LL  (\ 

-  (d  In  /d  In  4)  -  (d  In  t  /d  In  4)  (15) 

within  the  interval  Lq  <  L  <  L^.  For  x  and  x  exactly  (i.e., 
with  fixed  {3  and  y)  one  obtains  9  »  2(f3+y-2)  x  l/“  ^  ^  ^ “ ,  p  * 

((3-y-4)/(3-(3)  ,  and  v  -  ((3-3)/(p+y-2)  ,  in  which  case  f_(L)  is  given  exactly 
by  (13)  if  (3  >  3  and  (3  +  y  >  2.  The  subscripts  0  and  1  in  (13)  denote  evalu¬ 
ation  at  L  *  Lq  and  L  a  ,  respectively.  The  limit  r^(L)  —  a  in  (12)  yields 

f^(L)  ■  [ (4g-  4)/(4q“  C^)]?„(l^)  exactly.  Expansion  of  the  modified  Bessel 
functions  in  (13)  for  small  argument  yields  this  same  result  for  f„(L)  if  (3 
and  y  are  fixed  (i.e.,  independent  of  L) . 


The  above  development  of  a  modified  WKB  approximation  for  f  (L)  is  some- 

CO 

what  reminiscent  of  the  Green-^-Liouville  solutions  described  by  Jentsch  [1984], 
His  Green-Liouville  solutions  involved  hyperbolic-sine  (sinh)  functions, 
which  are  in  fact  proportional  to  modified  Bessel  functions  of  order  v  *  1/2. 
It  seems  that  (13)  is  the  appropriate  generalization  of  the  procedure  de¬ 
scribed  by  Jentsch  [1984]  to  situations  in  which  p  #  0. 


EIGENFUNCTIONS  AND  EIGENVALUES 


The  development  of  time-dependent  solutions  of  (1)  is  facilitated  by  the 
eigenfunction  expansion  shown  in  (11),  where 


L2(d/dL)[(D  /L2)(dg  /dL)]  -  r^g  (L)  +  \  g  (L) 

ll  n  q  n  n  n 


-1 


(d/d<t)[D$<5(dgn/d*)]  +  gn(L)  +  Vn(L)  =  °-  <16> 


The  eigenfunctions  gR(L)  are  required  to  vanish  both  at  L  =  LQ  and  at  L  *  L^. 


Eigenfunctions  corresponding  to  distinct  eigenvalues  X  and  X  are  necessari- 

n  m 


ly  orthogonal  in  the  sense  that 
L, 


f 


1  L"2  g  (L)  g  (L)  dL 
n  is 


nm 


0 ,  n  i  m 
1 ,  n  *  m 


(17) 


Given  a  complete  set  jg^tL)}  of  orthonormal  basis  functions  satisfying  (17) 
and  the  boundary  conditions  gn(Lp)  =  g^(L^)  *  0,  the  required  eigenfunctions 


gn(L)  and  eigenvalues  X^  can  be  obtained  by  diagonalizing  the  matrix  repre¬ 


sentation 


A 


nm 


£ 


L"2  i_  (L)  D  i  (L)  dL  + 
n  LL  m 


c 


L"2  g  (L)  *  ~l (L)  5  (L)  dL 

n  q  m 


(18) 


2  2 

of  the  transport  operator  A  =-  L  ( d/9L)  [  (D. .  /L  )(d/3L)]  +  (1/t  ).  An  optimal 


set  g^(L)  of  basis  functions  would  be  one  that  can  be  constructed  by  means 


of  a  fairly  simple  prescription,  but  one  that  makes  the  off-diagonal  elements 


of  A  especially  small  in  absolute  value, 
nm 

It  is  evident  from  (18)  that  all  the  eigenvalues  of  A  are  positive. 

nm 

This  situation  corresponds,  of  course,  to  temporal  decay  of  the  expansion 

coefficients  a  (t)  in  (II).  Moreover,  if  the  eigenvalues  \  are  ordered  (as 
n  n 

usual)  so  that  0  <  X.^  <  X.^  <  <  .  .  .,  then  it  follows  from  (18)  that  (with 

increasing  n)  radial  diffusion  becomes  increasingly  important  (compared  to 

charge  exchange)  for  the  determination  of  X.^.  The  presence  of  derivatives 

of  the  gn(L)  in  the  first  term  (but  not  in  the  second  term)  on  the  right-hand 

side  of  (18)  assures  this.  Except  for  the  term  t  ^g  (L),  which  has  no  coun- 

q  n 

terpart  in  their  paper,  the  eigenvalue  equation  specified  by  (16)  is  identi¬ 
cal  in  form  to  the  one  for  which  Schulz  and  Boucher  [1984]  successfully  con¬ 
structed  an  optimal  set  of  orthonormal  basis  functions  by  means  of  a  variant 
of  the  WKB  approximation.  A  further  variant  of  that  construction  is  required 
here,  since  the  boundary  conditions  of  the  present  radial-diffusion  problem 
differ  from  the  boundary  conditions  appropriate  to  the  pitch-angle  diffusion 
problem  treated  by  Schulz  and  Boucher  [1984]. 


and  (20)  is  given  by 


@  s 
n 


*  /yv-* 


where  9^  is  the  nCtl  positive  root  (n-0, l ,2, . . .)  of  C*(0^) 

C*(0°)  -  J  (0°)Y  (o0°)  -  J  (o0°)Y  (0°)  -  0, 

v  n  v  n  v  n  n  v  n 


(21) 

0,  i.e. ,  where 

(22) 


and  where 


*  0  *  0 

Since  C  (9  )  ■  C  (o0  )  »  0  for  each  value  of  n,  it  thus  follows  from  (19)- 
v  n  w  n 

(23)* that  gn(Lg)  *  8n(I^)  *  0»  as  is  required.  The  normalization  and  mutual 

orthogonality  [in  the  sense  of  (17)]  of  the  basis  functions  gn(L)  specified 

by  (19)  can  be  verified  by  using  9  as  the  variable  of  integration  and  in- 

n 

voicing  certain  indefinite  integrals  evaluated  by  Watson  [1944,  pp.  148-149], 
The  optimal  order  v  of  the  ordinary  Bessel  functions  in  (20a)  is  given  by  v  * 
l/(p+2),  where  p  is  a  representative  value  of 


p  s  (d  In  DLL/d  In  ;)  -  4(d  In  L/d  In  ()  *  d  In  D  /d  In  ;  (24) 

within  the  interval  L.  $  L  <  L.  (0^  >9  >  o0^).  For  D.  T  x  exactlv  (i.e., 

0  Inn  n  LL 

with  fixed  (3  >3)  one  obtains  0  *  (Ln/L)  ^  ^^0®,  a  *  (Ln/L.)^  p  * 

n  u  n  u  l 

(4-(3)/((3-3) ,  and  v  *  ((3-3)  /  (p— 2)  .  The  basis  functions  specified  by  (19) 

should  lead  to  a  nearly  diagonal  matrix  representation  A  ,  as  defined  by 

nra 

(18),  of  the  transport  operator  As-  (d/3<S)  [D^  0/9<X>)  ]  +  (1/t  ).  In  other 

words,  the  diagonal  element  A  (at  least  for  n  >  4)  should  greatlv  exceed 

nn 

the  absolute  value  of  each  off-diagonal  element  A  (=*  A  )  in  the  same  row 

6  nm  mn 

or  column  of  the  matrix.  This  major  benefit  of  the  WKB  construction  of  basis 

functions  enables  the  eigenvalues  and  eigenvectors  ot  A  (and  therefore  the 

a  nm 

eigenvalues  and  eigenfunctions  of  the  transport  operator  A  )  to  be  evaluated 
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by  means  of  a  rapidly  convergent  perturbation  theory.  The  formal  results 
[Schulz  and  Boucher,  1984]  are 


nn 


-  £ 


A  .  A. 
nk  kn 


(25a) 


k^n  ^kk  ^nn 


and 


*„(L)  '  Un„*n<W  +  E  Vl<l) ' 


(25b) 


kj*n 


where 


kn 


0  A  -  A.  . 

nn  nn  kk 


A.  + 
kn 


A..  A. 

_jy _ is. 


j^k.n  Ann  Ajj 


(25c) 


for  k  #  n  and 


U 


,  2, -1/2 


nn 


(25d) 


t1'  +  E  (Ukn/Unn>2'' 

ki*n 

to  assure  the  unitarity  of  the  transformation  from  the  jgn(L)|  to  the  jgn(L)}, 


i.e.,  to  assure  that  the  g^d)  are  likewise  normalized  in  accordance  with  (17). 


A  further  use  of  the  orthonormal  basis  functions  gn(L)  specified  by  (19) 


is  to  eliminate  altogether  the  presumably  small  discrepancy  between  the  exact 
steady-state  solution  f  (L)  of  (1)  and  the  modified  WKB-approximate  steady- 


state  solution  f  (L)  given  by  (13).  This  can  be  done  by  formally  expanding 


the  discrepancy  as  a  general  linear  superposition  of  the  g^(L)  and  inserting 


the  formal  expansion,  viz.. 


?  (L) 


f  (L) 


-  V  A  g 

nrm 


(L), 


(26) 


m 


into  (1)  for  3f/dt  »  0.  The  result  (after  the  usual  straightforward  steps) 


is  a  set  of  coupled  linear  equations  given  by 

L, 


Ya  a  -  f  L~2  i 

nm  m  J .  n 

_  “A 


(L)  Dtt  (df  / dL)  dL 

L.L  3 


L, 


-i, 


(L)  T  (L)  f  (L)  dL 

q  3 
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(27) 


for  the  expansion  coefficients  that  should  be  inserted  in  (26).  The  solu¬ 
tion  of  (27)  is  numerically  well-determined,  since  the  matrix  A  given  by 

nm 

(18)  is  supposed  to  be  nearly  diagonal  when  the  basis  functions  in(L)  are 
constructed  in  accordance  with  (19)-(24). 


APPLICATIONS 


This  report  was  previously  prepared  as  a  short  paper.  This  prevented 

the  inclusion  of  numerical  results  illustrating  the  usefulness  of  analytical 

methods  described  above  for  solving  problems  in  radial-diffusion  theory. 

However,  the  numerical  results  of  Schulz  and  Boucher  [1984],  showing  (for 

example)  that  the  off-diagonal  elements  of  A  are  consistently  smaller  (by 

run 

one  to  several  orders  of  magnitude)  in  absolute  value  than  the  corresponding 
diagonal  elements  when  analogous  methods  are  applied  to  a  prototypical  pitch- 
angle  diffusion  problem,  suggest  chat  the  present  approach  will  be  found 
highly  advantageous  when  applied  numerically  to  radial-diffusion  problems  as 
well.  Moreover,  the  ease  with  which  charge  exchange  can  be  incorporated,  as 
in  (19),  into  the  radial-diffusion  problem  suggests  that  charge  exchange 
could  similarly  be  handled  together  with  pitch-angle  diffusion  in  problems 
that  require  this,  e.g.,  in  studies  of  the  evolving  pitch-angle  distribution 
of  ring-current  ions  after  charge  exchange  has  made  the  equatorial  distribu¬ 
tion  anisotropic  enough  to  generate  electromagnetic  ion-cyclotron  waves  [cf. 
Cornwall,  1977].  For  this  latter  application  the  basis  functions  could  be 
constructed  according  to  the  prescription  of  Schulz  and  Boucher  [1984].  A 
long-range  goal  is  to  treat  the  simultaneous  occurrence  of  radial  diffusion 
and  pitch-angle  diffusion,  in  which  case  the  L-dependent  eigenvalues  of  the 
pitch-angle  diffusion  operator  will  presumably  enter  the  mathematical  de¬ 
scription  of  radial  transport  [cf.  Walt ,  1970]  in  somewhat  the  same  way  that 


the  charge-exchange  rate  r  1 (L)  enters  (1).  However,  the  bimodal  (radial/ 

’  12 


pitch-angle)  diffusion  problem  is  complicated  by  the  absence  of  a  kinematical 
quantity  that  both  modes  of  diffusion  simultaneously  conserve,  and  the  solu¬ 
tion  is  presumably  much  more  elusive  for  it  than  for  the  isolated  radial- 
diffusion  problem  treated  here. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  '’architect-engineer"  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exclmer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out -of -f ield-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  9olar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault -tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratoiy:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optic9 ,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  microwave/millimeter  wave  measurements, 
diagnostics  and  radiometry,  microwave/millimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric 
and  Ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effeccs  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


